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We investigate a scalar field dark energy model (i.e., 0CDM model) with massive neutrinos, where 
the scalar field possesses an inverse power-law potential, i.e., V(4>) oc 4>~ a ( a > 0). We find that 
the sum of neutrino masses Em„ has significant impacts on the CMB temperature power spectrum 
and on the matter power spectrum. In addition, the parameter a also has slight impacts on the 
spectra. A joint sample, including CMB data from Planck 2013 and WMAP9, galaxy clustering 
data from WiggleZ and BOSS DR11, and JLA compilation of Type la supernova observations, is 
adopted to confine the parameters. Within the context of the </>CDM model under consideration, 
the joint sample determines the cosmological parameters to high precision: the angular size of the 
sound horizon at recombination, the Thomson scattering optical depth due to reionization, the 
physical densities of baryons and cold dark matter, and the scalar spectral index are estimated to 
be 0* = (1.0415io;°on) x 1CT 2 , r = 0.0914^;^, D b h 2 = 0.0222 ± 0.0005, D c h 2 = 0.1177 ± 0.0036, 
and n 3 = 0.9644jtg'oii|, respectively, at 95% confidence level (CL). It turns out that a < 4.995 
at 95% CL for the <^>CDM model. And yet, the ACDM scenario corresponding to a = 0 is not 
ruled out at 95% CL. Moreover, we get Em„ < 0.262 eV at 95% CL for the <(>CDM model, while 
the corresponding one for the ACDM model is < 0.293 eV. The allowed scale of Em, in the 
0CDM model is a bit smaller than that in the ACDM model. It is consistent with the qualitative 
analysis, which reveals that the increases of a and Em„ both can result in the suppression of the 
matter power spectrum. As a consequence, when a is larger, in order to avoid suppressing the 
matter power spectrum too much, the value of E m v should be smaller. 

PACS numbers: 95.36.+X, 98.80.-k, 95.85.Ry, 13.35.Hb 


I. INTRODUCTION 

Neutrino is one of the important bonds linking nuclear physics, particle physics, astrophysics and cosmology [l[. In 
the Standard Model (SM) of particle physics, it is anticipated that there are three types, or “flavors”, of neutrinos: 
electron neutrino (p e ), muon neutrino ( v M ) and tau neutrino (i/ T ), which are also dubbed as three normal/active 
neutrinos. Besides that, neutrinos are assumed to be massless in the SM of particle physics 

It was first predicted by Bruno Pontecorvo in 1957 that if neutrinos are massive the neutrino flavor should be 
unstable, that is called neutrino (flavor) oscillations |3). Briefly put, neutrino oscillation is a phenomenon that a 
neutrino produced in a definite flavor is observed in a different flavor after traveling some distances. In other words, 
neutrinos are able to oscillate among the three available flavors while they propagate through space. Nowadays there 
are compelling evidences for neutrino oscillations from a variety of experimental data on solar, atmospheric, reactor 
and accelerator neutrinos. The discovery of neutrino oscillations implies that neutrinos have small but non-zero 
masses, with at least two species being non-relativistic today. However, the present experimental results on neutrino 
oscillations only measure the difference of two squared masses, such as = m 2 — m 2 and Am § 2 = m§ — m|, but 

give no hint on their absolute mass scales, mi, m 2 and m 3 are the neutrino mass eigenstates. For example, the solar 
neutrino analysis supplemented by KamLAND produces an estimate of Am^ ~ 8 x 10~ 5 eV 2 pjh, and the measurement 
of atmospheric neutrino oscillation by Super-Kamiokande I indicates A?n§ 2 ~ 3 x 10 _ 3 eV 2 Q. If it is the case of 
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oscillations among three light neutrinos, only two of the three Am| are independent, as + Am | 2 + Amf 3 = 0 , 

where A mf 3 = rri'j — m§. Recent reviews on progress in both theoretical and experimental aspects of neutrino 
oscillations can be found in [ 6 ]. 

A variety of cosmological tests are sensitive to the absolute scale of neutrino mass, such as the cosmic microwave 
background (CMB) radiation, galaxy surveys, and the Lyman-alpha forest (3]. In @, the effect of massive neutrinos 
on the Sunyaev-Zel’dovich and X-ray observables of galaxy clusters are investigated with a set of six very large 
cosmological simulations ( 8 h -3 Gpc 3 comoving volume). The analysis of current cosmological observations provides 
an upper bound on the total neutrino mass Y m„ (summed over the three neutrino families) of order 1 eV or 
less. However, the limits on Y m v from cosmology are rather model dependent and vary strongly with the data 
combination adopted. For example, in the framework of one-parameter extensions to the base ACDM model the 
Planck 2015 results @ give 95% upper limits on the sum of neutrino masses, i.e., Y m i / < 0-23 eV for a combination 
of Planck TT+lowP+lensing+ext, and Y m v < 0-59 eV for Planck TT, TE, EE+lowP+lensing, where “TT” denotes 
the combination of the TT likelihood at multipoles l > 30 and a low -l temperature-only likelihood, “TE” denotes the 
likelihood at l > 30 using TE spectra, and “EE” denotes the likelihood at l > 30 using EE spectra, “lowP” denotes the 
low -l Planck polarization data, “lensing” is the Planck lensing data, and “ext” represents the external data including 
the baryon acoustic oscillations (BAO), Type la supernovae (SNe la), and H 0 . In 0 , the power law and exponential 
types of viable f(R) theories along with massive neutrinos are studied. It shows that the allowed scales of Y m v 
in the viable f(R) models are greater than that in the ACDM model. The cases of fixing the effective number of 
neutrino species as N e g = 3.046 and treating N e g as a free parameters are both considered in 0. The former 
corresponds to just consider the active neutrinos without the effect of dark radiation. The latter corresponds to 
include the contribution of dark radiation (represented by AN e g = N e g — 3.046). For more details on dark radiation, 
we refer the reader to 0 - The model of holographic dark energy with massive neutrinos and/or dark radiation is 
investigated in jl2j , but the computed results from this model are not compared with those from the ACDM model. 
Actually, the ACDM model with massive neutrinos is discussed broadly with constraints from various cosmological 
observations fTdj . The time evolving of neutrino mass is also explored in the literature [l4| . For further details on 
neutrino cosmology, the reader is referred to recent reviews such as 0 , 0 - 

In this paper, we will discuss the constraints on the sum of neutrino masses Y m v i n the framework of </CDM model 
by using a combination of the CMB data from Planck 2013 and WMAP9, the galaxy clustering data from WiggleZ 
and BOSS surveys, and the JLA compilation of SNe la observations. The effect of dark radiation is not considered 
in this work, i.e., N e g = 3.046. We also assume that one of the three active neutrinos is massive, and the other two 
are massless. The </>CDM model — in which dark energy is modeled as a scalar field (j> with a gradually decreasing 
(in (/) potential V (</) — is a simple dynamical model with a slowly decreasing (in time) dark energy density. This 
model could resolve some of the puzzles of the ACDM model 0 ], such as the coincidence and fine-tuning problems. 
Here we focus on the scalar field with an inverse power-law potential V{(f >) oc <f>~ a , where a is a nonnegative constant 
00 - When a = 0 the </>CDM model is reduced to the corresponding ACDM case. The </>CDM model with this 
kind of V(<j) ) has been extensively investigated jTol - 10 . but without considering the massive neutrinos. 

The rest of the paper is organized as follows. In Sec.|TT]we present the background and perturbation evolutions of 
the </CDM model with massive neutrinos. The impacts of Y ^ and oc on the CMB temperature power spectrum 
and on the matter power spectrum are also discussed. Constraints from the cosmological data are derived in Sec. IIIII 
and the results for </CDM model are compared with those for the ACDM model. We summarize our main conclusions 
in Sec. [TV] 


II. THE 0CDM MODEL WITH MASSIVE NEUTRINOS 
A. Background evolution of the 0CDM model 


Quintessence as one of the popular scalar field dark energy models is a hypothetical form of dynamical dark energy 
to explain the late-time cosmic acceleration. Since quintessence is described by the scalar field </, the corresponding 
dark energy model can also be called as </>CDM model. In what follows, we will use the terms “quintessence” and 
“</>CDM” essentially interchangeably. We consider the self-interacting scalar field 4> minimally coupled to gravity on 
cosmological scales. The action of this </>CDM model is given by 


S = 


/ 




2 

— Ar . 

167T 


Crh + C 
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where g is the determinant of the metric g^, R is the Ricci scalar, m p = 1 /y/G is the Planck mass with G being the 
Newtonian constant of gravitation, C is the Lagrangian density for matter and radiation, and C^ is the Lagrangian 


3 


density for the field </>, given by 


C,h — 


p 

I67r 




( 2 ) 


where V(4>) is the field’s potential. In this work, we take a flat Friedmann-Lemaitre-Robertson-Walker (FLRW) metric 
for the background evolution, which is described by 

ds 2 = —dt 2 + a 2 5ijdx'dxd , (3) 

where x l is the conroving coordinate. a(t) is the scale factor usually normalized to unity now clq = a(z = 0) = 1 and 
related to the redshift 2 as a/ao = 1/(1 + z). Throughout, the subscript “0” denotes the value of a quantity today. 
By the variation of the action in Eq. © with respect to <j >, one can obtain the Klein-Gordon equation (equation of 
motion) for the scalar field 


(j> + 3 




( 4 ) 


For the </CDM model, there are many kinds of V(<f>) which can satisfy the requirement of the late-time accelerating 

In 1988, Peebles and Ratra [ 3 ] proposed a scalar field that is slowly rolling down 
at a large 0, where k and a are nonnegative parameters. This inverse power-law 


expansion of the universe 22]. 
with a potential V(<f>) = \nm 2 (j) 


2”""P 

potential can not only lead to the late-time acceleration of the universe but also partially solve the cosmological 
constant problems. The larger value of a induces the stronger time dependence of the scalar field energy density p^. 
When a = 0, this </>CDM model is reduced to the ACDM case. What is more, the parameter k depends on a (see 
EH HH for its dependence on a). 

The Friedmann equation of the 0CDM model with massive neutrinos can be written as 


H 2 { z ) — 7^2 (P b + Pc + P<j> + Pj + Pv), 


( 5 ) 


where pb, p c , Pc/> , p 7 and p v denote the energy densities of baryons, cold dark matter (CDM), scalar field dark energy, 
photons and neutrinos, and H(z) = a/a is the Hubble parameter. The energy density and pressure of the scalar field 
dark energy are given by 


9 

7TI 

P * = l6^ 2/2 + Vm 

and 

9 

777 ^ 

p * = 

Then, one can work out the equation of state (EoS) of the field <j>, 


— P(j>/p(j> — 


~ 2H(</>) 
</ 2 + 2 V (</>) ’ 


( 6 ) 

( 7 ) 

( 8 ) 


which is clearly bounded in the range — 1 < u)^ < 1 and usually non-constant. One can see that if the scalar field 4> 
rolls slowly enough such that the kinetic energy density is much less than the potential energy density, i.e. f (c/>), 
the pressure P of the field will become negative with —> —1. 

Based on Eqs. 0 and 0, along with the initial conditions described in Refs. 13 , one can numerically 
compute the Hubble parameter H(z). We also introduce the dimensionless density parameter for each component as, 
0.x = Px/Per, where the index “X” denotes the individual components, such as radiation (“r”), neutrino (V”) and 
matter (“m”). The critical energy density is expressed as p cr = 3H 2 m 2 /(8n). VL m is the energy density of matter 
including both baryons and CDM. fis the total neutrino energy density which scales as oc a~ 4 at early times, and 
thereafter evolves as oc a -3 after the non-relativistic transition. One can see that the massive neutrinos behave like 
the radiation at early times and like the matter later. 
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B. Cosmological perturbation of the 0CDM model 


Let us consider perturbations of the flat FLRW metric in the Newtonian Gauge 
perturbed metric is given by 


In this gauge, the linear 


ds 2 = a 2 (if) [—(1 + 24')d?7 2 + (1 + 2$)5ijdx l dx : ’] , (9) 

where the scalar perturbations are dominant over vector or tensor perturbations, and 77 = f a~ 1 dt is the conformal 
time. The Newtonian force 4' gives rise to the dynamics of the perturbed fluids, while the curvature perturbation 
measures the local energy density fluctuations. The linear perturbation theory is a good tool both for describing the 
early universe at any scales, and the recent universe on the largest scales. 

It has been shown in [25| that for self-interacting scalar field dark energy models it is phenomenologically sufficient 
to regard the dark energy component as a perfect fluid. We treat each component in the universe as perfect fluid, 
including the baryon, CDM, photon, neutrino and scalar field dark energy. In the perfect fluid approach, the perturbed 
Einstein equations lead to the following Eqs. (USD - m in the Fourier space: 

S' x + 3 Ha(cl x - uj x )6 x ,= -(1 + u x )(fi x + 3$') ( 10 ) 


Q\ 


x 


Ha( 1 — 3oj x ) + 


u x 


1 + U) X 


Ox = k 2 


i,X 


1 + LUX 


Sx + 4 > + o x 


( 11 ) 


= 47 xGa 2 p l [<5jc + 3 Ha(uj x + 1 )9 x /k 2 ] , (12) 

and 

4- = (13) 

The great advantage of linear theory is to obtain independent equations of evolution for each Fourier mode. All of 
the perturbed quantities (5 X ,9 X , 4/, d^efc.) are functions of space x and time t, where X denotes each perfect fluid 
composing the universe. In the linear perturbations, the anisotropic stress o x is negligible for the perfect fluids. 
Note that a prime represents a derivative with respect to the conformal time 77 . The spatial variation of density 
fluctuations is expressed by the density contrast S x = 5px/px = (px — Px)/px , and px is the background energy 
density of component X. In the approximation of negligible irrotational flow, the divergence of the peculiar velocity 
v x , 9x = V ■ v x can be used to describe the fluid motion. In the Fourier space , we have 0 x = ik ■ Vx- While 
u>x = Px/Px is the equation of state of each component, and c 2 x = 5P x /bpx represents the sound velocity. Eq. 
(flOll is called as the (perturbed) continuity equation , that states the conservation of local density. Eq. (fill) is called as 
the Euler equation, that represents the conservation of local energy momentum, and describes dynamics of perturbed 
fluids originated by the Newtonian force T. The curvature perturbation <f> is constrained to the local inhomogeneity 
via the Poisson equation Eq. m- We can get Eq. (THU) under the assumption that the perturbed fluid remains a 
perfect fluid. These equations (Hod - CSD completely determine the dynamical evolution of large scale structure (LSS) 
of the universe, within a given expansion history H. 


C. Matter power spectrum and CMB power spectrum in the 0CDM model 

In the framework of </>CDM model, we qualitatively investigate the impacts of parameters a and on the matter 
power spectrum and on the CMB power spectrum. The analyses are performed with the CAMB Boltzmann code f2(|. 

Neutrinos rarely interact with matter after thermal decoupling, so they are treated as free streaming particles. 
Massive neutrinos are the only particles that present the transition from radiation to matter. Before the non- 
relativistic transition the neutrinos behave like radiation. Thus, when the neutrino mass Sm„ increases, the time of 
radiation/matter equality is postponed gradually, and a eq increases. The value of T,m v can affect the matter power 
spectrum and the CMB power spectrum mainly resulting from a change in the time of equality, that provides a 
potential way to constrain it through CMB and LSS observations 7, 12, (27]. In Fig. [I] we show the impacts of 
neutrino mass Sm„ on the matter power spectrum P(k) and the CMB temperature anisotropy spectrum Cf T . The 
upper panels show the cases for <^>CDM model with varying values of Em„, where a is fixed as a = 1, and other 
parameters are fixed based on the recent Planck results jlj. For comparison, we also display the cases for ACDM 
model in the lower panels. For both <^CDM and ACDM models, the matter power spectrum is gradually suppressed 
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with the increase of Y>m u , however, the effect is more significant on small scales than that on large scales. One 
possible reason is that the neutrino perturbations do not contribute to gravitational clustering on scales smaller than 
the free-streaming scale, while on the very large scales neutrino perturbations are never affected by free streaming, 
and they become indistinguishable from CDM perturbations in the non-relativistic regime jj]. The CMB temperature 
anisotropy spectrum Cf T is insensitive to the variation of Y,m v in both ACDM and c/>CDM models. 

The parameter a indicates the dynamics of dark energy, and then it can affect the expansion history of the universe 
and the redshift of matter/dark energy equality. When a increases, the expansion of the universe occurs more rapidly, 
and the epoch of dark energy domination begins earlier [Tof . For these reasons, the variation of a can have signatures 
in the CMB map and the matter clustering. The impacts of a on P{k) and Cj T are presented in Fig. [2] We choose 
a = 0, 1, and 10 as examples, where a = 0 corresponds to the ACDM scenario. The values of other parameters are 
kept fixed. We find that P(k) is slightly suppressed with the increase of a, and the effect is a bit significant on large 
scales than that on small scales. The CMB temperature anisotropy spectrum C) TT is a little sensitive to the variation 
of a on the low-1 tail, which may arise from the late Integrated Sachs-Wolfe (ISW) effect. Anyhow, CMB and LSS 
observations are efficient to distinguish between ACDM and </>CDM models. 

III. OBSERVATIONAL CONSTRAINTS 

The observational data sets used to constrain the cosmological parameters are described as follows, including the 
galaxy clustering, CMB and SNe la measurements. 


A. Cosmological data sets 

1. Galaxy clustering measurements 

Galaxy clustering distilled from the galaxy redshift survey is powerful as cosmological probe |28j j, that can allow 
us to measure the cosmic expansion history through the measurement of BAO, and the growth history of cosmic 
large scale structure through measurements of redshift-space distortions. The length scale of BAO, the comoving 
sound horizon at the baryon drag epoch r s (zd ), can be applied as a cosmological ruler and accurately calibrated by 
observations of the CMB radiation. The position of the BAO peak in the angle-averaged galaxy clustering pattern is 
usually quantified in term of the volume averaged distance I29( 

D v =[(l + z) 2 D 2 A cz/H(z)] 1/3 . (14) 

It is common to report the BAO distance measurements as combinations of the angular diameter distance, Da{z), 
and the Hubble parameter, ff(z), such as 


^/ \ = H oV^ D v( z ) 

cz 

d z = D v (z)/r s {z d )- 


(15) 


(16) 


The redshifts of galaxies include indistinguishable contributions from both the Hubble recession and the peculiar 
velocity of the galaxies themselves, so that there are errors in the distances we assign to galaxies. The differences 
between the redshift-inferred distances and true distances are known as redshift-space distortions (RSD) [30|. In 
another word, the RSD are introduced in the observed clustering pattern by galaxy peculiar motions. As a consequence, 
the correlation function and the power spectrum measured in the redshift space are different from those in the real 
space, which have to be corrected to be expressed in real space. Because the effects of RSD couple the density and 
velocity fields, the RSD signals within the correlation function are difficult to model. On different scales, peculiar 
motions produce different types of distortions to the power spectrum. On small scales, i.e., in the cluster cores, the 
peculiar velocities of galaxies are almost randomly oriented, that cause the structures to a ppe ar elongated along the 
line of sight (LOS) when viewed in redshift space (i.e., the so called “finger of God” effect) |3lj], leading to a damping 
of the clustering. On large scales, because of gravitational growth, the galaxies tend to fall towards high-density 
regions, and flow away from low-density regions, such that the galaxy clustering in redshift space is enhanced in the 
LOS direction compared to the transverse direction [32|. The RSD effect on large scales can be described by linear 
theory (H, [H| , while the “finger of God” effect is a non-linear phenomenon. On large scales where the gravitational 
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growth is linear, measuring the relative clustering in both LOS and transverse directions leads to measurements of 
the parameter combination f (z e ff)as(z e g), where z e g is the effective redshift. / is the growth rate of cosmic structure, 
which is associated with the evolution of matter density perturbations S m via the relation / = dlnS m /d\na. In the 
linear regime, the linear growth rate can be expressed as / = dlnD(z)/dlna, where D(z ) = S m (z) / S m (z = 0) is 
the linear growth factor normalized such that D(z = 0) = 1. a${z) is the root-nrean-square amplitude of the matter 
fluctuations in spheres of 8 h~ x Mpc, and as{z = 0)/as{z) = D{z = 0 )/D(z). Thus, one can figure out 

f{z)a 8 (z) = erg (17) 

am a 

where Og = as{z = 0). In linear theory, the galaxy bias b and the growth rate / are degenerate with a%, so the 
RSD measurements are better presented in terms of b(z)a$(z) and f(z)a 3 (z), rather than /(z). Currently, the bias- 
independent parameter combination /(z)er 8 (z) measured by RSD are widely used [20l till l.'S5l. 

The Alcock-Paczynski (AP) test [36| is proved to be a significant link between BAO and RSD. AP test states that 
if an astrophysical structure is spherically symmetric or isotropic, then it should possess equal comoving transverse 
and radial sizes. An AP measurement is carried out by comparing the observed transverse and radial dimensions 
of objects. While the AP test is equally valid for an isotropic process such as the two-point statistics of galaxy 
clustering. The apparent anisotropy of the two-dimensional correlation function of galaxies mainly arises from the 
geometry and expansion of the universe which should be correctly embodied in the fiducial cosmological model, and 
the RSD effect which is supposed to be marginalized by using appropriate RSD model (see [30] for a recent review 
of RSD models). According to the requirement of AP test, the signature of BAO should have identical comoving 
sizes (i.e., r s ) in transverse and radial dimensions. The observed transverse dimension is the angular projection 
A 8 = r s /[(l + z)Da(z)]. The radial one is the redshift projection Az = r s H(z)/c. The relative radial/transverse 
distortion depends on the value of 


F(z) = Az/Ad 

= (1 + z)D A {z)H(z)/c, (18) 

where F(z) is dubbed as the AP distortion parameter. 

By combining the BAO peak, AP test and RSD effect, one can report the galaxy clustering effectively as joint 
measurements of (A, F, fas) or (d z , F, fas)- These joint measurements are extremely good at helping to constrain basic 
cosmological parameters and distinguish between the dark energy models. By using large-scale structure measurements 
from the WiggleZ Dark Energy Survey [37], Blake et al. (2012) [38| have performed joint constraints of (A, F, fas) in 
three overlapping redshift slices with effective redsliifts z e g = (0.44,0.6,0.73). Utilizing these data, it is straightforward 
to put constraints on the model parameters by calculating the corresponding XwiggieZ> given by 

Xwigglez = (*obs - X th )C~\X ohs - X th ) T (19) 

The observational data vector is 

A 0 bs = [Ai, A 2 , A 3 , Fx,F 2 , F 3 , fas,i, fas, 2 , fas, 3 ], (20) 

i.e., AAbs = [0.474,0.442, 0.424, 0.482,0.650,0.865,0.413,0.390,0.437] by using the maximum likelihood estimations of 
(A, F, fas) listed in Table 1 of (38|. The vector of theoretical values is 

X t h = [A(zi),A(z 2 ),A(z 3 ),F(z 1 ),F(z 2 ),F(z 3 ),fas(zi),fas(z 2 ),fas(z 3 )}, (21) 

where [zi,z 2: z 3 ] = [0.44,0.6,0.73], and the corresponding theoretical values of ( A,F,fa 8 ) can be obtained with 
Eqs. m, © and ED. respectively. C is a 9 x 9 covariance matrix between parameters and redshift slices, and the 
value of 10 3 C is listed in Table 2 of {3&J, that is achieved by generating 400 lognormal realizations for each WiggleZ 
survey region and redshift slice with the methods described in [34| . In the analysis of [38[, the fitting formulae provided 
by Jennings et al. (2011) (39] have been taken as the fiducial RSD model. The effect of different choices of the RSD 
model is also considered in Section 3.4 of [38|. It turns out that the systematic error induced from modeling RSD is 
much lower than the statistical error in the measurement. 

Joint measurements of (d z ,F, fas) at an effective redshift of z e s = 0.57 are provided in Samushia et al (2014) [40|] 
by utilizing the observed anisotropic clustering of galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS) 
Data Release 11 (DR11) CM ASS sample [ill] . We employ this data set in our analysis with the chi-squared statistic 

XbOSS = (bobs ~ bt h)Coil 1 (Fobs — l^h)^ 


(22) 
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The observational data vector is i^bs = [d z ,F,fcr$\, i.e., Fobs = [13.85,0.6725,0.4412] by using the mean values pre¬ 
sented in Eq. (30) of [40]. The vector of theoretical values is Y th = \d z (z e s), F(2 e ff), /er 8 (z ef f)], where the corresponding 
theoretical values of (d z ,F, fa 8 ) can be obtained with Eas. lfl^l) . (fISl) and m, respectively. The covariance matrix 
Cov of measurements is listed in Eq. (31) of (40[ . A suite of 600 PTHalo simulations are used to estimate the covari¬ 
ance matrix (see [42j for details of mock generation). In the analysis of [40|, the “streaming model”-based approach 
developed in [43| has been adopted to model the RSD signal, that has been demonstrated to fit the monopole and 
quadrupole of the galaxy correlation function with better than percent level precision to scales above 25/i~ 1 Mpc, for 
galaxies with bias of b cs 2. 

The galaxy clustering (GC) measurements from WiggleZ and BOSS DR11 are both employed in this study. Thus, 
the corresponding chi-squared statistic is expressed as 

XgC = XwiggleZ + Xboss> (23) 

where XwiggleZ an d Xboss are given by Eqs. (fl9l) and ([22]) . respectively. 


2. CMB power spectrum measurements 

The CMB radiation deemed as the afterglow of the big bang can supply us with some information of the very 
early universe. The observations of CMB provide another independent test for the existence of dark energy. The 
recent precise measurements of the CMB radiation from Planck and Wilkinson Microwave Anisotropy Probe (WMAP) 
projects can efficiently improve the accuracy of constraining the cosmological parameters. Currently, the Planck 2015 
results have come out but the Planck 2015 likelihoods are not yet available. Given this, we use the low multipoles 
(2 < l < 49) and high multipoles (50 < l < 2479) temperature power spectrum likelihoods from Planck 2013 fiij ]. 
together with the low multipoles (l < 23) polarization power spectrum likelihoods from nine-year WMAP (WMAP9) 
j45l |. To employ the previously mentioned CMB power spectrum data in the analysis, we compute the yjf MB statistic 

Xcmb = £(C? bB - Cf h )-M«7'(C? b8 - Cfl (24) 

ll' 

where C° hs is the observational value of the related power spectrum, Cj h is the corresponding theoretical value in 
the framework of the cosmological model under consideration, and Ai is the covariance matrix for the best-fit data 
spectrum. 


3. Magnitude-redshift measurements of Type la supernovae 


The first direct evidence for the cosmic acceleration came from SNe la observations, which provide the measurement 
of the cosmic expansion history through the measured luminosity distance as a function of redshift, di(z) = (1 +z)r(z). 
In the spatially flat universe, the comoving distance r(z) from the observer to redshift z is given by 


r(z; p) 


c A dz' 
H~oJo E(z'-,p)’ 


(25) 


wherein p denotes the parameter space of the considered cosmological model, and E = H/Hq is the dimensionless 
Hubble parameter. 

Here, we use the “joint light-curve analysis” (JLA) compilation of SNe la [46j, which is a joint analysis of SNe la 
observations including several low-redshift samples (z < 0.1), all three seasons from the SDSS-II (0.05 < z < 0.4), 
three years from SNLS (0.2 < 2 < 1), and 14 very high redshift (0.7 < 2 < 1.4) from the Hubble Space Telescope 
(HST) observations. It totals 740 spectroscopically confirmed SNe la with high quality light curves. In Ref [461 ]. 
SALT2 light-curve model [13, EH have been used to fit the supernova light curves of the JLA sample. From the 
observational point of view, the distance modulus of a SN la can be yielded from its light curve with an empirical 
linear relation: 


M ° bs =m* B - {M b -axli +/3 x C) (26) 

The light-curve parameters (nig, Xi,C) result from the fit of the SALT2 light-curve model to the photometric data, 
where m* B corresponds to the observed peak magnitude in rest-frame B band, X\ describes the time stretching of the 
light-curve, and C describes the supernova color at maximum brightness, a , f) and M B are nuisance parameters in 



the distance estimate, which are estimated simultaneously with the cosmological parameters and then marginalized 
over when obtaining the parameters of interest, wherein Mg is the absolute B-band magnitude. The theoretical 
(predicted) distance modulus is 

P,Mo) = 51og 10 [D L (2;p)] +/i 0 , (27) 

where [ 1 q = 42.38 — 51og 10 h , which is also treated as a nuisance parameter, and the Hubble-free luminosity distance 
is given by 

dl{z • p > s v <ii(2) = (i+z) L ifc) ■ (28) 

The best-fit cosmological parameters from SNe la data are determined by minimizing 
740 

XSNe = Pi M b) - M th ’ l (^5 P, Mo) 

i,j =1 

where Cov is the covariance matrix of data vector /7^ s . The values of the covariance matrix Cov and the SALT2 fit 
parameters {m* B , Xi,C) are available from Ref. [461] . 


Cov. 


-l 


^(a, 13, Mb) - V th ' 3 fa; p,Mo) 


(29) 


B. Results and analysis 

In our analysis, the likelihood is assumed to be Gaussian, thus we have the total likelihood 

£ oc e _Xtot ^ 2 , (30) 


where Xtot is constructed as 

Xtot — XGC + XCMB + XSNei (31) 

wherein Xq C , Xcmb an d XsNe are gi yen by Eqs. (ESI- (EH and (OH1) . respectively, and denote the contributions from 
galaxy clustering, CMB and SNe la data sets described above. 

We derive the posterior probability distributions of mrameters with Markov Chain Monte Carlo (MCMC) explo¬ 
ration using the February 2015 version of CosmoMC [4{j. The parameter space of the ACDM model is 

Pacdm = {n b h 2 ,n c h 2 , 100 9mc, T, ln(10 10 A s ), n s , Em,}, (32) 

where Q b h 2 and £l c h 2 , respectively, stand for the baryon and CDM densities today, Bmc is an approximation to 
d* = r s {z£)/D a{z*) (i.e., the angular size of the sound horizon at the time of decoupling) that is used in CosmoMC 
and is based on fitting formulae given in [50j . r refers to the Thomson scattering optical depth due to reionization, 
n s and A s are the power-law index and amplitude of the power-law scalar primordial power spectrum of curvature 
perturbations, and E m u is the sum of neutrino masses. The parameter space of the <)>CDM model is 

P^cdm = {STfo/z, 2 , S7 c /i 2 , 1OO0 mc> T , ln(10 10 A s ), n s , E m v , a}, (33) 

which has one more parameter than that of ACDM model, where a determines the steepness of the scalar field 
potential in the framework of </>CDM model. 

The one-dimensional (ID) probability distributions and two-dimensional (2D) contours for the cosmological param¬ 
eters of interest are shown in Fig. [3] for ACDM model and in Fig. [3] for 0CDM model. It shows that constraints from 
the joint sample are quite restrictive, though there are degeneracies between some parameters, such as the degen¬ 
eracies in the fl m — Hq and erg — Yiniv planes. In addition, the differences between the marginalized likelihoods and 
the mean likelihoods are modest in ID and 2D plots. It implies that the distributions of the parameters are almost 
Gaussian. We also present best-fit values and mean values with 95% confidence limits for the parameters of interest 
in Table [T] both for ACDM and <jiCDM models. We find a < 4.995 at 95% CL for the </>CDM model, while the ACDM 
scenario corresponding to a = 0 is not ruled out at this confidence level. The constraints on Vl b h 2 , Q. c h 2 , 1OO0mc> D 
ln(10 lo A s ), n s , D m , er 8 and Hq are consistent at 95% CL for these two models. 

Here, we pay attention to the constraints on T,m u . Note that Em y is in unit of eV. The best-fit vale is E?n„ = 
0.038(0.043) with Em„ < 0.262(0.293) at 95% CL in the framework of 0CDM (ACDM) model. The allowed neutrino 
mass scale in the <^CDM model is a bit smaller than that in the ACDM model. As it is shown in Figs. [T] and [2] the 
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increases of a and T,m u both can result in the suppression of the matter power spectrum. Therefore, when a is larger, 
in order to avoid suppressing the matter power spectrum too much, the value of E m v should be smaller. Consequently, 
the case with a > 0, i.e., 0CDM model, has smaller E m v \ correspondingly, the case with a = 0, i.e., ACDM scenario, 
has larger T,m„. Additionally, in Ref. 0, they obtained Em„ < 0.200 eV at 95% CL in the ACDM model, that is 
consistent with our result. With the results presented in Table II of 0 , one can see that our constraints on f Ibh 1 2 3 4 5 6 7 , 
D c ft 2 , t and n s are also consistent with theirs at 95% CL for the ACDM model. 


IV. CONCLUSION 

We have concentrated on a quintessence model (or called as 0CDM model) of dark energy with massive neutrinos. 
In the </>CDM model under consideration, the scalar field (f> is taken as a candidate of dark energy to drive the late¬ 
time acceleration of the universe with an inverse power-law potential V(<p) oc </>““ (a > 0). The larger value of a 
corresponds to the stronger time dependence of the scalar field energy density. When a = 0, it is reduced to the 
corresponding ACDM scenario. The linear perturbation theory is employed in the framework of this model. Through 
qualitative analyses, we find that the increases of the sum of neutrino masses Em„ and the parameter a both can 
gradually suppress the matter power spectrum P(k). It implies that when the value of a is bigger, in order to avoid 
suppressing the matter power spectrum too much, Etoj, should be smaller. It is in accordance with the results from the 
observational constraints. The variations of these two parameters also can have signatures in the CMB temperature 
anisotropy spectrum Cj T . In order to make a comparison, the impacts of E m v on P(k) and Cj T in the context of 
ACDM model have also been presented. 

A combination of the CMB data from Planck 2013 and WMAP9, the galaxy clustering data from WiggleZ and 
BOSS DR11, and the JLA compilation of the SNe la observations is used to constrain the parameters. The results 
indicate that constraints on the cosmological parameters from this joint sample are quite restrictive. It turns out that 
Em„ < 0.262 eV (95% CL) in the framework of 0CDM model and Em„ < 0.293 eV (95% CL) in the ACDM model. 
The allowed neutrino mass scale in the c/>CDM model is a little shrunk comparing to that in the ACDM model. In 
Ref. 0, it is concluded that the allowed neutrino mass scales in the viable f(R) models are bigger than that in the 
ACDM model. Given this, we can infer that the allowed scale of E m v in our d>CDM model must be smaller than those 
in the viable f(R) models. In addition, we get a < 4.995 at 95% CL for the <(>CDM model, meanwhile, the ACDM 
scenario corresponding to a = 0 is not ruled out. Consequently, the observational data that we have employed here 
still cannot distinguish whether dark energy is a time-independent cosmological constant or a time-varying dynamical 
component. 
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FIG. 1: Impacts of the sum of neutrino masses E rth, on the matter power spectrum P(k ) and on the CMB temperature power 
spectrum Cf T in the 0CDM (upper panels) and ACDM (lower panels) models. E m v is varied, and other parameters are kept 
fixed. 
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FIG. 2: Impacts of the parameter a on the matter power spectrum P(k) and on the CMB temperature power spectrum Cj T 
in the framework of 0CDM model, a is varied, and other parameters are kept fixed. 


Parameters 
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0CDM Model 

Best-fit values 
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TABLE I: Fitting results from the joint sample. We present the best-fit values (i.e., the parameters that maximize the overall 
likelihood), and the mean values with 95% confidence limits for the parameters of interest. Where Em„ is in unit of eV, and 
Ho is in unit of km/s/Mpc. The top block contains parameters with uniform priors that are varied in the MCMC chains. The 
lower block defines various derived parameters. 
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ACDM Model 





FIG. 3: The ID and 2D probability distributions of parameters of interest in the ACDM model constrained with the joint 
sample. In the ID plots, the solid lines denote the marginalized likelihoods and the dotted lines correspond to the mean 
likelihoods. In the 2D plots, the contours refer to the marginalized likelihoods while the colors refer to the mean likelihoods. 
The contours correspond to 68%, 95% and 99% confidence levels. 
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FIG. 4: The ID marginalized distribution and 2D contours of parameters of interest in the c/>CDM model constrained from 
the joint sample. The implications of line styles and colors are the same as those in Fig. [3j 







































